International Journal of Power Electronics and Drive Systems (IJPEDS) 
Vol. 14, No. 4, December 2023, pp. 1965~1975 
ISSN: 2088-8694, DOI: 10.1159 1/ijpeds.v 14.14.pp 1965-1975 O 1965 


Evaluation of the possibility of chaos for doubly-fed induction 
generator in wind power generation system 


Co Nhu Van!”, Nguyen Phung Quang’, Nguyen Thanh Hai! 
‘Faculty of Electrical-Electronic Engineering, University of Transport and Communications, Hanoi, Vietnam 
"Institute for Control Engineering and Automation, Hanoi University of Science and Technology, Hanoi, Vietnam 


Article Info ABSTRACT 

Article history: The doubly-fed induction generator (DFIG) has the advantages of 
; compactness and low cost, so it has been applied more and more in practice. 

Received Dec 14, 2022 However, the structure of the DFIG system is complex, multivariable, 

Revised Jun 11, 2023 nonlinear, and difficult to control. Chaos is a phenomenon that only occurs 

Accepted Jun 25, 2023 with nonlinear dynamical systems, sensitive to initial conditions and aperiodic 


but it is governed by deterministic laws, unlike random perturbations. For 
most powertrain systems, chaos causes disadvantages, reduce performance, 
Keywords: and even destruction of the system. By method of analysis based on theory 
and simulation, the paper will demonstrate the chaotic phenomenon for DFIG 


Bifurcation when the parameters of the system change into a certain value. From that set 

Control chaos . of parameters, the controller is designed to avoid or eliminate that chaotic 

Doubly fed induction generator behavior, so that the system works more reliably and with better quality. 

Lyapunov exponents 

Portrait phase This is an open access article under the CC BY-SA license. 
©) BY SA 

Corresponding Author: 

Nguyen Thanh Hai 


Faculty of Electrical-Electronic Engineering, University of Transport and Communications 
No. 3 Cau Giay Street, Ha Noi, Vietnam 
Email: nguyenthanhhai@utc.edu.vn 


NOMENCLATURE 
R,,R; : Rotor, stator resistance isastsqs tras trg :dq components of the stator, rotor 
current 
Lm, L:, Ls : Mutual, rotor, stator inductance Us, Uy : Vector of stator, rotor voltage 
T;, Ts : Rotor, stator time constant Urq) Urq» Usa» Us: dq components of the stator, rotor 
voltage 
Co) : Total leakage factor W,, ¥, : Vector of rotor, stator flux 
W, W;,@s : Mechanical rotor velocity, rotor and Y.q, Yq : dq components of the stator flux 
stator circuit velocity 
ig, ly : Vector of stator, rotor current Pra, Pra : dq components of the rotor flux 


1. INTRODUCTION 

DFIG in the power generation system has the outstanding characteristics that the stator side is directly 
connected to the grid and the rotor side is connected to two back-to-back voltage source converters to the grid. 
Because the control device is located in the rotor, the power of the control device is much smaller than the 
generator power (in the limited speed range, the power of the converter is only 30% of the power transmitted 
to the grid [1], [2]), which is economically attractive, especially when the generator power is large. The system 
is capable of operating in a slip coefficient in a fairly wide range (+33% compared with the synchronous 
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speed) [3], [4] allowing to make good use of the power source hybridized by the main machine, which is 
working in either over synchronous or sub synchronous mode. In both modes, the machine supplies active 
power P to the grid on the stator side, and on the rotor side, the machine draws reactive power Q from the grid 
in the sub synchronous mode and returns the reactive power Q to the grid in the over synchronous mode. 

Due to the great benefits of DFIG, this configuration has become now very popular for variable-speed 
wind turbines, so there are many control methods for DFIG to perfect the system, typically studied [5]—[10]. 
However, in industrial practice, reliability is extremely important, therefore, it is necessary to study more 
deeply and widely to help the system achieve the best efficiency and quality. The parameters of DFIG can 
change with temperature, life, and load conditions. Due to the challenging working environment of wind farms. 
DFIG is vulnerable to faults like gearbox faults, power converter faults, stator winding faults, rotor winding 
faults and velocity sensor faults. The word [11]-[13], from which the system may fall into a state of chaotic 
work, which results in poor system working quality and is the reason for issues and failures. 

Recently, the research investigation into the bifurcations and chaos of motor drives and permanent 
magnet synchronous generators has received much attention [14]-[19], enabling a better comprehension of the 
system's dynamical activities, which may provide some information that is helpful for actual control and design 
work. However, research on chaos for DFIG is only a few scattered studies, that do not cover all issues of chaos 
in DFIG. While DFIG is considered complex, the parameters of DFIG can vary with temperature, life, 
environment, load conditions, etc., so the system may fall into a chaotic behavior, which makes the system's 
working quality poor, which is the main cause of system failure. Thus, it is necessary to have an overview and 
study more fully, applying chaos to control DFIG for wind power generators will lay the foundation for new 
research, thereby providing appropriate control structures to complete the system. 

The paper is organized as follows: i) The overview of the chaos phenomenon and the identification to 
control chaos will be presented in section 2; ii) Section 3 will present some mathematical models of DFIG as 
a basis for implementing section 4 — chaos in DFIG, to perfect the system, a controller is introduced to eliminate 
or avoid the phenomenon of chaos; iii) The discussion will be done in section 5; and iv) Finally, section 6 
presents the conclusions the of study. 


2. CHAOS AND APPLICATIONS IN CONTROL 

In the real world, any system is disturbed by external noise, the oscillations of systems are always 
nonlinear, the linearization of oscillation processes is nothing more than a simplification of the real oscillation 
problem, and in computer simulations, a small perturbation appears due to numerical round-off. A little 
perturbation that occurs in a chaotic system accumulates exponentially over time, dramatically altering the 
system's behavior [20], [21]. There are many chaotic regimes in nature and in human-made devices, chaotic 
dynamics is one of the most general methods of the evolution of nonlinear systems. When chaos intensifies 
mixing and speeds up chemical reactions, it is advantageous because it creates a strong mechanism for the 
transfer of mass and heat. Chaos is a common undesired phenomenon that can, for instance, lead to additional 
mechanical fatigue of the elements of construction due to their irregular vibrations. In a chaotic regime, there 
is a chance that no resonant energy absorption will push the system parameters over safe limits. 

Chaos brings undesired results and can be self-maintaining, affecting the quality of the system. This 
phenomenon only manifests in nonlinear systems, is nonperiodic, sensitive to initial conditions, and follows 
certain laws that are not the same as noise. It was not until the last decades of the 20th century that chaos theory 
began to be deeply explored for powertrain systems. Because the input and output variables of the powertrain 
vary from the time, most of the calculations are instantaneous, to simplify the control problem we have to 
accept the idealized input variable and treat the drivetrain parameter as constant. However, that turns them into 
noise, which causes chaos in the system for a while. With the control requirements set out in the drive array, 
the accuracy of response, and stable operation, the application of chaos theory about the control object to 
increase reliability is very feasible. 

The most crucial aspect of chaos motion is that it is highly sensitive to the initial conditions of the 
system. That is, a very small difference in the input that is amplified will make a very large difference in the 
output. This is shown very clearly by the simple system of (1) of Lorenz-1963 [22] when the initial conditions 
are different (Figures 1). 

The parameters have the values o = 24.1, p = 10, B = 8/3, Figure 1 shows that, with very little different 
initial conditions, from the value [0 1.79980000000000000002 0] to the value [0 1.7998000000000001 0], but 
the state of the system has changed from a steady state (Figure 1(a)) to a chaotic state (Figure 1(b)). 
Basically, a chaotic system has the following properties [23]-[25]: 

— Nonlinearity: Chaos only occurs in nonlinear dynamical systems or systems that exhibit a certain degree of 
nonlinearity. 
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— Determinism: chaos can be predicted by simple deterministic equations, determining the domain of 
parameters corresponding to the chaotic solution. That is, chaos obeys one or more deterministic equations, 
without an element of randomness, or probability. 

— Sensitive to initial conditions: A small change in the system's initial state can result in drastically different 
final-state behavior. As a result, even though the behavior of the system is determined by deterministic 
fundamental laws, it is impossible to anticipate it over the long term. 

— Aperiodicity: Chaotic orbits are aperiodic, but follow a certain law or principle. People are often interested 
in a special form of this family called the attractor. The orbits wander forever in an internal domain, without 
reaching a fixed point or a closed orbit, they are attracted to a complex geometrical body called the strange 
attractive, which is a chaotic phenomenon [20]. In addition to being a strange mesomorphic structure, the 
chaotic attractor is also dynamically strange. The orbits are attracted to the attractive, but they are unstable 
on that set and are sensitive to initial conditions — a characteristic of chaos motion. Once inside the strange 
attractor, the phase portraits are confined to it and can approach any point arbitrarily attractor to the suction 
set but never repeat the same at a later. The phase trajectory is not stable anywhere on the strange attractor, 
but overall, the attractor is very stable as shown in Figure 2. 


ao) 

dy _ = = 

o = x(p 2) -y () 
dz 

rien aa Bz 


Figure 2. The attractor from the simulation of Lorenz equations on MATLAB with o = 28, p = 10, B = 8/3 
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— Defining Chaos: We must first recognize that chaos never exists in a linear dynamic system. Thus, when 
we refer to chaos, we actually mean nonlinear systems. Chaotic motion is not present in all nonlinear 
systems, though. The Fourier analysis method, time responses, phase portraits, bifurcation diagrams of the 
maximum value of the state variable over time, and calculations of Lyapunov exponents are a few of the 
techniques used to identify chaos. These techniques are used to show the chaotic behavior when changing 
the characteristic value, to be able to identify the chaotic working area of the object, that propose controls 
method to eliminate the self-sustaining oscillations with high amplitude and irregular variation, return the 
system to a stable working state. As analyzed, the DFIG system has a complex structure, especially DFIG 
in wind powers working in unstable conditions from wind energy, which easily leads to a chaotic state. 
Therefore, identification and elimination from the "third state" in wind power generation systems using 
DFIG to increase reliability are important goals of this research. 


3. STATE EQUATION OF DFIG 

Decoupled regulation of active and reactive current components is one of the key control goals. This 
advises using a reference frame for the control architecture that is stator voltage orientated. The real axis of the 
grid voltage vector un might be chosen as the d axis in generating systems like wind power plants where the 
stator is directly connected to the grid. The following relations of the DFIM are obtained [26]: 


Ls » . 
Wa = i, tsa + tg = 0 
1 _ ds, oe GG (2) 
Wsq — Ln 154 au lrg a IWs| 


where ¥! = Y/Lm. 


tid = een 2) ete (2 fa — Oy) + Fo ttng = 54 
<a — (2 + ~*) ing = Orlra + = (= Vea + wa) + stra -_ sq (3) 


Where o = 1 - I2,/(L,L,), T; = L,/Rs, T, = L,/R,.In grid voltage orientated coordinates: u 
0; Usa=Us; Pog = 0; Yo = YY and: 


sq 


sin @ = ae M¢g = —=2 a Weqird (4) 

The rotor current component i,g acts as a control variable for the generator torque mG and 
respectively for the active power P, because the stator flux Y; is governed by the grid voltage and can be seen 
as constant. The power factor cos@ or the reactive power Q can be controlled by the control variable i,,. The 
system of (3) shows very clearly the multivariable and nonlinear of the system, which is a necessary condition 
for the system to be chaotic. Thus, when the DFIG system operates under certain conditions, it can lead to 
chaotic behavior. This issue will be considered in detail in content 4 and content 5 below. 


4. CHAOTICS PHENOMENA IN DFIG 
4.1. Some causes of DFIG working unstable, bifurcation, and chaotic behavior 
— When the controller parameter changes 

From the classical control structure PI for DFIG in [27], based on theoretical analysis, from the 
Jacobian matrix, it is possible to determine the conjugate complex solutions whose real parts correspond from 
negative to less negative and become positive, as shown in the Table 1. One pair of complex eigenvalues' real 
parts become less negative as Kp,, increases, and at a critical value of Kp,,, the real parts switch from negative 
to positive, leading to Hopf bifurcation. 
— When the mechanical torque changes 

The research of Yang et al. [28] has determined that the DFIG wind generator is capable of bifurcation 
to instability when the mechanical torque Ti, increases to characteristic values through simulation and 
theoretical analysis. On the basis of theoretical analysis, the loci of the Jacobian's eigenvalues are computed 
and the analysis shows that the system loses stability via a Hopf bifurcation, there exist four pairs of complex 
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conjugating eigenvalues as T,, is varied. From Table 2 we see, for a small T,,, value, all four pairs have a negative 
real part; as Tm increased, the real part of the eigenvalues pairs became less negative and at the critical value 
of T,, the real part changes from negative to 0; As T,, continued to increase, the real part became positive and 
the system became unstable. This suggests that the instability is caused by Hopf bifurcation. 
— When DFIG is connected to a “weak” grid with unbalanced loads and parasitic 

At low wind speeds, bifurcation has been identified in the DFIG wind turbine, which manifests as a 
low-frequency oscillation of the DC link capacitor voltage of the system when the DFIG wind turbine is 
coupled to a "weak" grid that has parasitic and unbalanced loads, and if the wind speed further increases, 
chaotic behavior may emerge [29]. 
— When the parameters of DFIG vary 

This is acommon phenomenon for DFIG during operation, and there have been studies demonstrating 
that chaos occurs when the parameters change to a certain value [30]—[32]. In [31], from system (5). 


= Tha + Gig 

4 34 Big +B 

dé q awe (5) 
ee 

aE = O(ig — @ 


Where ig and i, are the transformed direct-axis and quadrature axis stator current respectively; tg, ig are the 


qd 
transformed direct-axis and the quadrature-axis stator voltage components respectively; @ is the transformed 


angle speed of the motor; o and y are system parameters. In the case of ig(0) = &,(0) = T,(0) =0; o = 5.46; 
Y = 20; (ig, ig, ©)" = (0.01, 0.01, 0.01)" the system is chaotic. 


Table 1. The eigenvalues for increasing value of Kp, (T;,, = 0:00076) [27] 
Kpo Eigenvalues Comment 
128 -1.0626 + j292.65, -33.743 + j281.95, -10.638 + j25.341, -3.0069 + j24.433 Stable 
129 -0.2951 + j292.74, -34.711 + j282.86, -10.638 + j25.341, -3.0069 + j24.433 Stable 
129.4 0.00136 + j292.78, -35.087 + j283.21, -10.638 + j25.341, -3.0069 + j24.433 Unstable 
130 0.43483 + j292.87, -35.641 + j283.72, -10.639 + j25.341, -3.0069 + j24.433 Unstable 


Table 2. Eigenvalues for increasing value of T,, [28] 


Tn Eigenvalues Comment 
-0.3  -0.0654+ j0.956, -0.0492+ j0.875, -0.0315+ j0.110, -0.0115+ j0.0851 Stable 
-0.4 = -0.0341+ j0.931, -0.0785+ j0.902, -0.334+ j0.127, -0.0116+ j0.0851 Stable 
-0.5 -0.0018+ j0.934, -0.1096+ j0.902, -0.0347+ j0.145, -0.0116+ j0.0851 Stable 

-0.507 0+ j0.935, -0.111+ j0.902, -0.0348+ j0.153, -0.0116+ j0.0851 Hopf bifurcation 
-0.6 _0.022+ j0.938, -0.1325+ j0.903, -0.0355+ j0.250, -0.0117+ j0.0852 Unstable 


The most recent study on chaotic phenomena and chaos control for DFIG is found in [32] typical 
research work, in the normal mode, through theoretical analysis of 0-1 test the system was determined to be 
stable. When the winding fault occurs in the system, the resistance value decreases accordingly Rs = Rr = 
0.02 2 (in nominal mode Rs = Rr = 1.2) At the same time, when the rotor speed sensor fails wo = 2 w, now the 
phenomenon of chaos has appeared. A fault observer based on the DFIG model is constructed to monitor the 
output current and the speed of the DFIG system in order to better watch the output of DFIG. Accordingly, the 
new sliding mode control's design is simpler, and it allows for better observation of the control effect. However, 
the condition of simultaneous failure of the rotor winding, stator winding, and speed sensor is rare in practice. 
To clarify the phenomenon of chaos for DFIG and methods to control the chaos. The following will do through 
the rotor winding fault condition. 


4.2. Demonstration of DFIG's chaotic phenomenon based on theoretical analysis and simulation 
4.2.1. State equation of DFIG 
From the equation of motion of the rotor [27]. 


“= "(7,-1,-2w) (6) 


dt J . 


where: T, = 5 Mphm Gaia - Uni lea)e The control method is based on the grid voltage orientated coordinates, 


we have Yq =0, Yq = Y;, From (2) we can rewrite in (6) as (7). 
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dw _ Np (3 Nplm¥s . D 
Go ee oe (7) 
at J \2 Ls Np 


From the first 2 equations of the system of (3), combined with (7) we get the system of equations 
representing the rotor current on the rotating dq coordinate and the rotor speed of DFIG as (8): 


di 1/1 1-o\. 3 1-o Wy, 1 1-o 

rd = _ ( + tra ial ti Mia wo +—u,g - Us 

dt o \Tr Ts o Lm oLr oLm 
dirg 1 (2 1-0) . ; 1-o 1 1 

= —-(-+—")i,, -(@, - wig + — = +—u 

dt o\T, | Ts J 74 (ws Nira G Tslm | oly ™4 (8) 
dw Np (3NpLmYs . D 
— plat af (ail Aan vd —- —w - T, 
dt J \2~ Ls np 


where: Yq = 0, Yq = Y, Usq = Us =U, Usg =0, Wig = Yq/Ln, Psa = Vsa/Lm, @s= 2mf. The (8) is 
reduced to (9): 


Xy = Cy xX1 + (Ws — Xz) XQ — CoX3 + CzUrg — C4 
Xq = CyX2 — (Ws — X3)Xy + C5 + CE6Urg (9) 


Xz = C7Xy — CgXz3 — Col) 
é ‘ 1/1 1-0 1-o 1 1-o 1-o 1 
where: Xi =1,93 X2=1,93X3=W3C1=— ( + )se2= W303 =—3 Ca = u,3 C5 = W.5c6 =—. 
1 rd» A2 rq? *3 > 1 o \Tp Ts >» &2 oLim go & Ly’ 4 oLm go &S oLmTs so 6 Ly 
Ny 3NpLmvs  3nZLmY D n 
p 3 NplLm?s pemts, . D.. _™p 
C7 = = 3;Cg=—;Co= 
J 2 Ls 2JLs J J 


4.2.2. With the normal operating mode of the system 

Table 3 lists the key system parameters that we employed in our simulations. Figure 3 and Figure 4 
illustrate when the system is operating under normal conditions. Figure 3(a) is time series diagrams of i,q, 
Figure 3(b) is phase trajectory diagram between i,q and w. Figure 4(a) is phase trajectory diagram i,q, tg and 
w@, Figure 4(b) illustrate all 3 Lyapunov exponents have values are negative. 
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Figure 3. The system is stable under normal conditions (a) time series diagrams of i,q and 
(b) phase trajectory diagram between i,q and w 
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Figure 4. The system is stable under normal conditions (a) phase trajectory diagram i,q, tg and w and 
(b) variation of Lyapunov's exponent 


Table 3. Main parameters of the DFIG system 


Symbol Parameters Value Symbol Parameters Value 
U Grid voltage (V) 575 TL load torque (N.m) 0 
f Grid frequency (Hz) 50 Rg, Ls Resistance (2) and inductance (uH) of stator 20, 68 
J Moment of inertia (kgm2) 2 R,, L; Resistance (42) and inductance (uH) of rotor 21, 69 
Np Pole pairs 3 Lin Mutual inductance (wH) 68 


D Damping coefficient (Nm/rad/s) 0 


4.2.3. The chaotic phenomenon is determined based on theoretical analysis 

In the event of a rotor winding failure, which results in the rotor resistance and inductance changing 
to a different value (other system parameters remain unchanged): R; = 15 0; L, = 38.6 wH. The equilibrium 
points of the system are determined by solving the system of (9), getting the following results: 


= _ Cg 
x= aa 

7 

= Co Cs 
X_ = (W, — X3) —T, - = 
2 ( Ss ae L Cy (10) 

CoTL — Cs 2Ws5CoTL\ — CoTL W5Cs5 CzCoTL 
Ott z+ (S- g— SS *) e+ e— =F - o,+ =0 
C407 C407 C407 Cy C7 


Substituting the parameters in (10) and solving, we get an equilibrium point: E(O; -16310; 314,2). 
Stability evaluation of equilibria through Jacobian matrix: From (9) we get: 


Cy (w, — X3) —C2—X, 
Je = |—(@s — *3) Cy Xy 
C7 0 0 
A- Cy —(@,— %3) (2 +%x, 
det|Al — J; | =|(w, — %3) A-G —X, 
—C7 0 a 


=A - 2c,A? + [cf + (Ws — 3)? + C7(Cg + X2)]A— C107 (Cg + Xz) - C7%1 (Ws — %3) 


The eigenvalues of the Jacobian matrix calculated at the equilibrium point are what determine the 
stability of the equilibrium point. It is customary to solve equation det[AI-J;] = 0 for A. 


A3 - 2c, A? + [c?2 + (W, — %3)? + C7(Cp + JI A- C1 C7 (Cp + X>)- C7X,(W, -— %3)=0 (11) 
With the equilibrium point E(0; -16310; 314,2) the solution of (11) is determined: 
A, = 1,57 + j716,84; A, = 1,57 - 7716,84; Az = -1,04 


We see, A,,A2 are conjugate complex solutions with positive real part, so the equilibrium point E is 
unstable. Thus, the system has 1 unstable equilibrium point. Next, we determine the Lyapunov exponent to 
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determine whether the system is chaotic or not. From the system of (9) with the given set of parameters, we 
can determine that the system has 3 Lyapunov's exponent (A,,A,A3) with the same value as Table 4. We see, 
the system always exists at least | Lyapunov's exponent that is always positive, so the system has a chaotic 
phenomenon. For more clarity, the following simulations are performed in the following content. 


4.2.4. The chaotic phenomenon is determined through simulation 

From the set of parameters of the system given in section 4.2.3. The simulation results show that the 
rotor speed fluctuates violently and does not follow any rules as shown in Figure 5(a), component amperage i,q, 
i-g oscillates around the equilibrium point, the trajectory does not repeat the past trajectory as shown in Figure 5(b) 
and Figure 6(a). At the same time there is always a positive Lyapunov exponent as shown in Figure 6(b). 


Table 4. Lyapunov's exponent 


Time A Ay a3 Time AL Ap 13 
0.1 -6.624482 4.647723. 4.647723 0.7 6.133618 -1.428603 —_-2.596103 
0.2 2.705223 —--17.230366 +: 16.653823 0.8 9.357681 —--1.139491 6.118786 
0.3 3.217109 -3.772286 ~—-2.672108 0.9 7.246686 ~—--1.128662_—_-4.022807 
0.4 2.930204 —--2.546148 ~—1.732098 1.0 9.052879 —_-0.836013—_-6.123998 
0.5 8.012185 —-1.881995 —_-4.016063 1.1 9.919040 -0.756869 —--7.072680 


0.6 1.044099 -0.624790 1.691378 
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Figure 5. The system works chaotically when an error occurs (a) time series diagrams of w and (b) phase 
trajectory diagram between i,q and w 
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Figure 6. The system works chaotically when an error occurs (a) phase trajectory diagram between i,q, irq 
and w and (b) variation of Lyapunov's exponent 
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5. ELIMINATE CHAOTIC PHENOMENA TO STABILIZE THE SYSTEM 

The chaotic phenomenon for DFIG is when the parameters of the system change to a certain value, 
which has been clearly shown through simulation and theoretical analysis (as shown in section 4.2). Chaos is 
harmful to the DFIG system, so it is necessary to design the controller to avoid or eliminate it when it occurs. 
The content below will take care of those issues. 

When chaos occurs, the controller's task is to bring the trajectory to the equilibrium point, we use the 
following control law: 


Uy = —ky(%, — X1); Ug = —k2 (x2 — X2); ug = —k3 (x3 — X3) 

Ky = CX, + (Wy — X3)X2 — C2X3 — C3 — ky (X41 — X41) 

X_ = —(Ws — X3)X1 + Cy X2 + Cy — kp (x2 — X2) 

X3 = C7X1 — Col, — k3(X3 — X3) (12) 


where X1, X2, X3 are equilibrium points determined from (10). 


A-Cy +k, —(w, =o X3) C2 + Xz 
det|Al —Je| =| (@,- %3) A-cq tk, —-x 
—Cz 0 At+k3 


=/A3 + (—2c, + ky + ky + kz)d? + [(W, — %3) 7% +7 (Cp + XQ) + kg (—2c, + ky +k) A+ 
[—cyc7(cp + %2) — C7X (Ws — ¥3) + c7kg (co + Xp) + k3(ws — X3) 7] (13) 


According to the Hurwitz stability criterion, from the conditions for the system to be stable, we can 
choose the parameters: k, = 520,000; kz = 1; k3 = 1:1) Control to eliminate the chaotic phenomenon. When 
chaos occurs, at 0.4 seconds the controller is inserted. Through simulation, we can see that the controller has 
eliminated the chaotic phenomenon, bringing the system to stable operation as shown in Figure 7(a); and 
ii) Control to avoid chaos: Once the boundary has been determined, the DFIG system can fall into a chaotic 
working state. The chaos avoidance controller will be included in the system from the very beginning. Based 
on the controller has been defined and the parameters have been selected. The simulation results show that the 
system operates stably as shown in Figure 7(b). 
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Figure 7. Simulation results when the controller is inserted for (a) time series diagrams of w and 
(b) time series diagrams of i, 


6. CONCLUSION 

Due to DFIG's many advantages, variable-speed wind turbines are now increasingly common in this 
design. However, research on chaos for DFIG is only a few scattered studies, that do not cover all issues of 
chaos in DFIG. therefore, it is necessary to study more deeply and widely to make the system more complete. 
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The paper has summarized the errors that can cause DFIG to fall into the phenomenon of bifurcation 
and chaos, to get a general picture. At the same time, through analysis on a theoretical basis and through 
simulation, which has demonstrated the chaos for the system when the rotor winding fails, the working area 
causing chaos of the system can be determined, thereby facilitating system parameterization to build a 
controller to avoid or eliminate that chaos, making the system work more reliably and with better quality. 
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